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Effects of interplanetary shock inclinations on nightside auroral power intensity 
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We derive fast forward interplanetary (IP) shock speeds and impact angles to study the geoeffec¬ 
tiveness of 461 IP shocks that occurred from January 1995 to December 2013 using ACE and Wind 
spacecraft data. The geomagnetic activity is inferred from the SuperMAG project data. SuperMAG 
is a large chain which employs more than 300 ground stations to compute enhanced versions of the 
traditional geomagnetic indices. The SuperMAG auroral electroject SME index, an enhanced ver¬ 
sion of the traditional AE index, is used as an auroral power (AP) indicator. AP intensity jumps 
triggered by shock impacts are correlated with both shock speed and impact angle. It is found 
that high AP intensity events typically occur when high speed IP shocks impact the Earth’s mag¬ 
netosphere with the shock normal almost parallel to the Sun-Earth line. This result suggests that 
symmetric and strong magnetospheric compression leads to favorable conditions for intense auroral 
power release, as shown previously by simulations and observations. Some potential mechanisms 
will be discussed. 

Keywords: Space physics, ionosphere-magnetosphere interaction, plasma physics. 


I. INTRODUCTION 

Interplanetary (IP) shocks result from the interaction 
of solar disturbances with the ambient solar wind [4p| . 
As they propagate throughout the heliosphere, IP shocks 
eventually interact with different bodies in the solar sys¬ 
tem, such as planets, moons, and even asteroid5 s. In the 
eventual cases in which IP shocks strike Earth, they inte¬ 
ract with the Earth’s magnetosphere, causing disturban¬ 
ces that can be detected in the near-Eartli space environ¬ 
ment, the whole magnetosphere and even the ionosphere. 
The first dramatic magnetospheric effect associated with 
IP shock impacts is the SSC/SI + (storm sudden com¬ 
mencement/sudden impulse), resulting from the sudden 
magnetospheric/magnetotail compression and the Earth¬ 
ward motion of the Chapman- Ferraro current. Other 
effects following shock-related SSC/SI+ events may also 
occur: geomagnetic storms (l3| . radiation belt perturba¬ 
tions [la,[53: and GICs (ground-induced currents). GICs 
may pose risks to electric power transmission systems le¬ 
ading to power grid disruptions and serious economic los¬ 
ses [ 3 , [lOi LLZt Sa] ■ GICs are also associated with corrosion 
of pipelines and their control systems [Hj].IP shocks are 
well known to sometimes trigger substorms as well fiH ]. 

IP shocks are well known to sometimes t rigg er subs- 
torrns as well d ! I [DIM EMI EE MM- In the 
early days before shock detection in interplanetary space, 
SSCs/SI + s were used to imply the impingement of an IP 
shock or tangential discontinuity (TD) onto the magne¬ 
tosphere as an attempt to explain geomagnetic activity 
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following SSC/SI + events. As shown later by [jjj, who 
statistically examined SSC/SI + events, intense au- ro- 
ral activity always occurred when SSC/SI + amplitudes 
were greater than 40 nT. EH later showed that most of 
SSCs/SI + s were caused by IP shocks rather than TDs. 
Precursor IMF B z events ~1.5 hr prior to shock arrival 
have bee used to identify when shocks would be geoeffec- 
tive and when they would not be jE 2j]. ;5l|, [54 - t5fil | 

Another important factor of IP shock geoeffectiveness 
is the IP shock impact angle, which is the angle between 
the shock normal vector and the Sun-Earth line. IP 
shocks with shock normals almost aligned with the Sun- 
Earth line are typically driven by CMEs (coronal mass 
ejections) EE whereas inclined shocks tend to be driven 
by CIRs (corotating interaction regions) more frequently 
[38l E3 • Several studies addressing geomagnetic activity 
following inclined IP shocks have been done in the past. 
For example, a longer than usual SSC rise time caused 
by the impact of an IP shock with the Earth’s magne¬ 
tosphere was observed by Wind spacecraft EE- These 
authors argued that this effect should be related to the 
high inclination of the shock normal in the equatorial 
plane. Due to its high impact angle, the IP shock took 
a longer time to sweep over the magnetosphere, com¬ 
pressing it gradually leading to a slow magnetospheric 
response. The authors suggested that simulation studies 
with the impact of inclined IP shocks on the Earth’s mag¬ 
netosphere should be carried out. This suggestion was 
taken by other authors [15,] who simulated the impact of 
two similar IP shocks, one inclined and the other frontal, 
on the Earth’s magnetosphere. They found that both 
systems evolved to very similar final quasi-steady states, 
although in the inclined case the system took more time 
to evolve in relation to the head-on case. Very similar 
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results were found with MHD numerical simulations by 
other studies as well [M M, Ml- A statistical study 
with more than 300 fast forward IP shocks was perfor¬ 
med [ 53 } to study the same effects, i.e., the effect of IP 
shock inclinations on the SSC rise time followed by IP 
shock impacts. They found that IP shocks impacting 
the Earth’s magnetosphere with shock normals almost 
aligned with the Sun-Earth line caused short SSC rise 
times when compared to inclined IP shock events. Their 
highest correlation occurred when the shocks were strong 
and almost head-on. These simulation and statistical re¬ 
sults confirmed the observation and suggestion made pre¬ 
viously by observations [49[ and numerical simulations as 
well [Td SUM HI- 

Impact of IP shocks on the Earth’s magnetosphere 
were investigated in other simulation studies SJ S3 with 
an MHD code especially developed to study effects gene¬ 
rated on the magnetosheath S3- In this case, the shock 
impact was frontal and the magnetospheric response was 
symmetric. Later, using the same MHD code, one of 
the previous authors a simulated case of the interaction 
of a similar IP shock, but with an inclined shock nor¬ 
mal, with the Earth’s magnetosphere. He found that the 
inclined shock took more time to travel through the mag¬ 
netosheath in comparison to his previous studies. This 
author argued that inclined shocks with large downs¬ 
tream v y component may generate asymmetries on both 
dawn and dusk sides leading to non-symmetric magne¬ 
tospheric compressions. Such compressions can lead to 
different SSC amplitudes depending upon the side of im¬ 
pact. He then suggested that such effects should be de¬ 
tectable by magnetometers on the ground. 

More recently, the impact of IP shocks with different 
shock normal inclinations on the Earth’s magnetosphere 
were simulated [3P, 33j. Using the OpenGGCM MHD 
code jUj], these authors simulated three different cases, 
namely, two inclined and one frontal. The second in¬ 
clined shock was twice as stronger as the other inclined 
shock, and the frontal shock had the same strength as 
the first inclined shock. The shock normals of the two 
inclined shocks lay in the meridian plane. The same 
authors reported that the head-on shock was more geo¬ 
effective than the inclined shocks, even more geoeffective 
than the strong inclined shock. For example, the frontal 
shock triggered substorm signatures and high nightside 
auroral energy dissipations not seen in the inclined ca¬ 
ses. They suggested that the frontal shock, whose shock 
normal was aligned with the Sun-Earth line, compressed 
the magnetotail symmetrically on both north and south 
sides. Such condition created an ideal scenario for the 
energy stored in the magnetotail to be leaked away and 
trigger auroral substorms at Earth [M Ml- These re¬ 
sults suggested the same authors to look for these effects 
in satellite and geomagnetic data. 

The sequence of this simulation work was conducted 
by the same authors {35j, j36j, who performed a statistical 
study of IP shock properties at 1 AU using a shock list 
with events from January 1995 to December 2013. They 


found that the yearly number of IP shocks is well correla¬ 
ted with solar activity, confirming previous observations 
Although shocks occur more frequently during so¬ 
lar maxima due to the higher occurrence of CMEs, CIRs 
tend to drive most shocks in solar minima. Due to the 
fact that CMEs tend to drive frontal shocks and CIRs 
tend to drive more inclined shocks, more frontal events 
tend to occur during solar maxima. They also reported 
that the majority of shocks found in the heliosphere at 
1 AU are weak shocks, with Mach numbers less than 3. 
In the same statistical study, they studied the effects of 
IP shock impact angles on substorm strength indicated 
by an enhanced version of the AL index. They found 
that almost frontal shocks were generally more geoeffec¬ 
tive than inclined shocks. Their strongest correlation was 
found in the case when shocks with high speed, or strong 
shocks, impacted the Earth almost frontally. 

The goal of this paper is to study geomagnetic acti¬ 
vity triggered by IP shocks in correlation to IP shock 
speeds and impact angles. Here the geomagnetic activity 
is represented by auroral power (AP) intensity as infer¬ 
red from an enhanced version of the auroral electroject 
index AE. Correlations are obtained from shock speed 
and impact angles using the same IP shock list published 
by a shock statistical study (3j|], the most extensive fast 
forward shock study done to date. In the following, in 
section 2, we present the data. In section 3, we report 
our results, which are summarize and briefly discussed in 
section 4. 


II. DATA ANALYSES 

Our data analysis is based on a list of 461 fast forward 
IP shocks from Wind and AGE data from January 1995 
up to December 2013 published by a recent statistical 
study [13 ■ The shock normal orientations were obtai¬ 
ned from different shock normal determination methods, 
such as the well known magnetic and velocity coplana¬ 
rity methods [5(, and the formulas that mix plasma and 
IMF data [l,, ,2{. All shocks were req uired to satisfy the 
Rankine-Hugoniot conditions [Hj, l50| . 

We use the SuperMAG geomagnetic station data to 
identify auroral power associated with shock impinge¬ 
ment. SuperMAG [HI is an international collabora¬ 
tion with a chain of more than 300 ground stations 
used to compute the SME, SMU, and SML indices 

f .28j, the enhanced versions of AE, AU, and AL 
respectively. The SuperMAG indices are very si¬ 
milar to the traditional I AG A indices: their main dif¬ 
ference is the fact that the former are computed based 
on data of a larger number of ground stations in com¬ 
parison to the latter. The SuperMAG data were obtai¬ 
ned from the websites http://supermag.jhuapl.edu/ and 
http://supermag.uib.no/. Technical issues related to Su¬ 
perMAG data analysis and assimilation were detailed in 
a paper entirely devoted to this subject [l3|. 

The SME index is used as a proxy for aurora power 



3 


(AP) determinations. This choice was based on a rela¬ 
tion found by other authors |28|. These authors calibra¬ 
ted the SME index with both Polar UVI instantaneous 
images and DMSP instantaneous maps to obtain possible 
correlations between SME and AP. Due to time resoluti¬ 
ons issues, the most relevant correlation found by them 
was between SME and AP as determined by Polar UVI. 
The linear relationship found in this previous work and 
used here is: 

AP = 0.048 x SME + 0.241 x (SME) 1 / 2 , (1) 

where AP is represented in GW, and the square root por¬ 
tion comes from the monoenergetic auroral contribution. 
In equation (1), AP was integrated over the northern 
hemisphere polar cap between 1800-0600 magnetic local 
time and 60° and 80° magnetic latitude. More specifi¬ 
cally, expression (1) indicates the nightside AP intensity 
as calculated from the SuperMAG SME index. Later, 
the SME index was confirmed to be the best choice to 
predict AP intensity instead of SMU and SML [29|. 



UT 


Figura 1: ACE observation of an interplanetary shock on 18 
April 2001 at 0005 UT and its consequent geomagnetic acti¬ 
vity. Top panel shows increase in dynamic pressure pv 2 (DP, 
in nPa) observed by ACE. Nearly 50 minutes later, as shown 
in the middle panel, the Earth’s magnetopause is struck by 
the IP shock, as indicated by SuperMAG ground stations with 
the sharp increase in the SMR index, in nT. Finally, the bot¬ 
tom panel shows the increase in AP(GW) approximately 35 
minutes after shock impact. 


The methodology used to record geomagnetic activity 
followed by IP shock impacts is shown by Figure 1. This 
figure represents an IP shock from the shock list compi¬ 
led by a previous statistical study [Hj. At 0005 UT on 
18 April 2001, ACE observed a sharp jump in the dy¬ 
namic pressure DP = pv 2 upstream of the Earth. After 
approximately 50 minutes, the IP shock impinged on the 
Earth’s magnetopause, and a sharp jump in SMR, the 
SuperMAG measurement associated with the ring cur¬ 
rent [26j, was recorded by SuperMAG ground stations. 
Then, in the next 40 minutes, the SuperMAG ground 
stations registered a peak in the SME index, which was 
used to plot AP in the last panel of Figure 1 according 
to equation (1). For all events in our statical analysis, 
the maximum measurements in AP followed by IP shock 
impacts were recorded in the time lag of 2 hours after 
shock impacts. If there are more than one AP peak in 
this time interval, the first one is chosen as the maximum 
associated with the IP shock. More details can be found 
in a previous work [35). 

III. RESULTS 

Figure 2 represents the statistical results obtained from 
our 461 IP shock events. Figure 2(a) shows the distri¬ 
bution of 9 Xn , the shock impact angle between the shock 
normal and the GSE Sun-Earth line. Angles close to 180° 
indicate that the shock normal vector is almost parallel to 
the Sun-Earth line. Most shocks had 9 Xn , indicating that 
they typically range from moderately inclined to almost 
frontal shocks. Figure 2(b) represents the distribution of 
shock speed vs in our shock list. The average shock speed 
is about 500 km/s in the Earth’s frame of reference, and 
most shocks have vs below the average. This result indi¬ 
cates that the IP shocks observed in the heliosphere at 1 
AU are predominantly weak IP shocks with Mach num¬ 
bers between 1 and 3. More details about the statistical 
results of the IP shocks in our database can be found in 
a previous paper [35 . 


Correlations of variations in auroral power, AAP, in 
GW, with the shock speed v s , in km/s, is shown in Fi¬ 
gure 3. For this parameter selection, the impact angle 
9 Xn is held in constant intervals while the shock speed is 
allowed to vary. The data are binned in three different 
categories: Figure 3(a), 120° < 9 Xn < 140°, highly incli¬ 
ned shocks; 2(b), 140° < 9 Xn < 160° , moderately incli¬ 
ned shocks; and 2(c), 160° < 9 Xn < 180°, almost frontal 
shocks. Here we consider events with low auroral activity 
when AAP < 20 GW, and events with high auroral acti¬ 
vity when AAP l 80 GW. Events with moderate auroral 
activity are between these two limits. Figure 3(a) shows 
that most highly inclined shock events with low auroral 
activity are associated with weak, low speed ( v s < 450 
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Figura 2: Statistical results of fast forward IP shocks observed 
by Wind and ACE from January 1995 to December 2013. 
Figure 1(a) shows the shock impact angle distribution, and 
Figure 1(b) represents shock speed (km/s) distribution. 


km/s) shocks. Strong, high speed (v a > 550 km/s) shocks 
are related to events with moderate auroral activity, with 
only one event that has low auroral activity being cau¬ 
sed by a strong shock. Events with moderate auroral 
activity are associated with all shock strength categories 
with approximately the same likelihood. There are no 
events with high auroral activity triggered by highly in¬ 
clined shocks in our database. The correlation coefficient 
in this case is R = 0.45 and the average of AP is 25.98 
GW. 

The intermediate category of shock strength has the 
largest number of events, as seen in Figure 3(b). In this 
case, all events with low auroral activity are triggered by 
weak or low speed shocks. Most events with moderate ac¬ 
tivity are associated with weak or moderate shocks. All 
events with high auroral activity are triggered by high 
speed shocks. The correlation coefficient is R = 0.55 and 
AP = 41.37 GW. Figure 3(c) shows that all weak auro¬ 
ral activity events (only three cases) are related to weak 
shocks. Events with moderate auroral activity are mos¬ 
tly associated with weak or moderate shocks, but some 
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Figura 3: Auroral power amplitude as a function of shock 
speed binned in three different shock impact angle categories: 
(a), 120° < 9 Xn < 140°; (b) 140° < 0 Xn < 160°; and (c), 160° 
< 9 X „ < 180° . 


are related to strong shocks. All events with intense au¬ 
roral activity are triggered by either moderate or strong 
shocks. The correlation coefficient R = 0.70 and the ave¬ 
rage AP = 64.09 GW are the highest in this category. 
These results are summarized in Table 1. 


A comparison among the three cases described above 
shows that, on aver- age, AAP increases with shock speed 
when the impact angle is close to 180o.The correlation co¬ 
efficient between the shock speed and AAP also increases 
for almost frontal shocks. 
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The opposite analysis is made in Figure 4, i.e., where 
the shock impact angles are allowed to vary keeping the 
shock speed binned in constant intervals. The three ca¬ 
tegories are: Figure 4(a), 300 < v s < 450 km/s, weak 
shocks; 450 < v s < 550 km/s, moderate shocks; and v s 
l 550 km/s, strong shocks. Figure 4(a) shows that weak 
shocks are associated with events with either weak or mo¬ 
derate auroral activity, and are not related to events with 
intense auroral activity. There are only a few weak highly 
inclined shocks, and most of them cause events with mo¬ 
derate auroral activity. Only a few strong highly inclined 
shocks cause events with low auroral activity. The cor¬ 
relation coefficient for highly inclined shocks, R = 0.39, 
and the average of 26.55 GW, are the lowest in this case. 
In the category of moderate shocks, the correlation is 
stronger, with R = 0.48, and the average is higher, with 
AP = 46.56 GW. There are only a few events with low 
auroral activity, and most of them are triggered by highly 
inclined shocks and just a few by inclined shocks. Mode¬ 
rate almost frontal shocks, seen in Figure 4(b), triggered 
either moderate or strong auroral activity events. Typi¬ 
cally, events with moderate auroral activity are triggered 
by moderate, strong, and weak shocks. There are only 
a few events with high auroral activity, and all of them 
are triggered by moderate almost frontal shocks. Finally, 
correlations for strong shock are represented by Figure 
4(c). Generally, strong shocks do not cause events with 
low auroral activity, with an exception of only one event 
caused by a highly inclined shock. Events with moderate 
AP activity are typically caused by inclined shocks, but 
they can also be triggered by highly inclined or almost 
frontal shocks. Events with intense auroral activity 235 
are caused mostly by almost frontal shocks, but a few 
events are caused by inclined shocks. The correlation co¬ 
efficient and AP average for strong shocks are the highest 
in this category, R = 0.79 and AP = 62.88 GW. Table 
1 summarizes the results obtained for correlations with 
shocks in all categories. 
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The analysis of the three panels in Figure 4 leads to 
a similar conclusion obtained in the analysis of Figure 3: 
strong, high speed shocks are generally much more geoef¬ 
fective than weak slow speed shocks, and their geoeffecti¬ 
veness increases if the IP shock impacts more frontally on 
the Earth’s magnetosphere. These general results were 
predicted previously in global MHD simulations [33| and 
confirmed experimentally with spacecraft and geomagne¬ 
tic data observations 1351. 


IV. SUMMARY AND CONCLUSIONS 

We have studied 461 fast forward interplanetary (IP) 
shocks using Wind and ACE satellite data from January 
1995 to December 2013. We correlated IP shock impact 
angles with geomagnetic activity (auroral power inten¬ 
sity) triggered by IP shock impacts. The primary result 


Figura 4: Auroral power amplitude as a function of shock 
impact angle binned in three different shock speed categories: 
(a), 300 < v s < 450 km/s; (b) 450 < v 3 < 550 km/s; and (c), 
v s > 550 km/s. 


Fixed impact angle (9 In , changed shock speed v 3 

category highly inclined moderately inclined almost frontal 


R 

AP 

0.45 

25.98 

0.55 

41.37 

0.70 

64.09 

Fixed shock speed v s , changed impact angle 9 Xn 

category 

weak 

moderate 

strong 

R 

0.39 

0.48 

0.79 

AP 

26.55 

46.56 

62.88 


Tabela I: Summary of the results obtained for the shock speed, 
shock impact angle, and AAP correlation analyses. 
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obtained here was that high speed shocks with shock nor¬ 
mal aligned along the Sun-Earth line (head-on shocks) 
cause the greatest auroral power release. The correlation 
coefficient for the cross correlation analysis in this case 
was 0.79, the highest of any performed in this study. This 
result confirms previous numerical simulation results [33}, 
in which frontal shocks led to stronger geomagnetic ac¬ 
tivity in comparison to the cases of inclined IP shocks. 
Observational results were reported in a subsequent work 
[35j], whose authors performed a statistical analysis cor¬ 
relating substorm strength and IP shock impact angles. 
In the case of fast (strong) shocks, events with the stron¬ 
gest geomagnetic activity occurred in the cases in which 
the shocks impacted the magnetopause almost frontally. 

To explain the above results, it should be first noted 
that shock compression of the magnetosphere is most 
effective when the inclination angle is frontal. Both 
the magnetosphere and magnetotail will be compres¬ 
sed the most for this orientation. Greater tail lobe fi¬ 
elds will require stronger cross tail currents to maintain 
them. Magnetosphere/magnetotail compression will lead 
to more flattened tail closed field lines. Slrock-triggering- 
substorm mechanisms were previously discussed by seve¬ 
ral other authors [j, [ 5 J (5J, [53] ■ Both current dis rup tion 
f22[ [23l . [ 37 I ] and magnetic reconnection pjJ [231425} are 
viable under these above conditions. 

The present results indicate the role of shock speed 
and inclination angle in geoeffectiveness of magnetosphe- 
ric energy release (auroral power). Thus this is another 
factor besides magnetospheric priming that must be ta¬ 
ken into account in assessing auroral power release. 
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